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On the Propagation of Waves through a Stratified Medium, with 
Special Reference to the Question of Reflection. 

By Lord Eayleigh, O.M., F.E.S. 
(Keceived December 8, 1911,— Eead January 11, 1912.) 

The medium is supposed to be such that its properties are everywhere 
a function of but one co-ordinate x, being of one uniform quality where x is 
less than a certain value x 1} and of another uniform quality (in general, 
different from the first) where x exceeds a greater value x m -i ; and the 
principal problem is the investigation of the reflection which in general 
ensues when plane waves in the first medium are incident upon the strati- 
fications. For the present we suppose the quality to be uniform through 
strata of finite thickness, the first transition occurring when x = x±, the 
second at x = x 2 , and the last at x = x m -i. 

The expressions for the waves in the various media in order may be taken 
to be 

<£i = Aie i ^ ct+h ^~ a ^ x ~ x ^^l\e i ^ ct+b ^ +a ^ x ~ x ^ > (1) 

d>2 = A2e i t ct+b y- a z( x ~ x M + 'B 2 e i [ ct+b y +a 2( x -zd] ) 

d>S = A3 e i ^ ct + h y~ a z( x ~ %>3 -f. g 3 e i [ct + hij + a 3 (x- 3:2)] 

and so on, the A J s and B J s denoting arbitrary constants. The first terms 
represent the waves travelling in the positive direction, the second those 
travelling in the negative direction ; and our principal aim is the determina- 
tion of the ratio B1/A1 imposed by the conditions of the problem, including 
the requirement that in the final medium there shall be no negative wave. 

As in the simple transition from one uniform medium to another (' Theory 
of Sound/ § 270), the symbols c and b are common to all the media, the first 
depending merely upon the periodicity, while the constancy of the second is 
required in order that the traces of the various waves on the surfaces of 
transition should move together — equivalent to the ordinary law of refrac- 
tion. In the usual optical notation, if V be the velocity of propagation and 
6 the angle of incidence, 

c = 2 ttV/ X, I = (2 tt/X) sin 0, a = (2 tt/\) cos 0, (2) 

where V/X, X" 1 sin 6 are the same in all the strata. On the other hand a is 

variable and is connected with the direction of propagation within the 
stratum by the relation 

a — b cot 6. (3) 
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The as are thus known in terms of the original angle of incidence and of 
the various refractive indices. 

Since the factor & ( ct+b ri runs through all our expressions, we may regard 
it as understood and write simply 

fa = Aie-^fr-^ + Bi^*-*!), (4) 

Hn. ... 

fa = A 2 e- ia ^ x ~~^ + B 2 e ia ^ x - x J, (5) 

fa = A s e~- ia ^ x -~^ + ^B se ^(x-x 2 ) f ^ 

<£ m = Awr^^-^i) + B^^-^-i). (7) 

In the problem of reflection we are to make B m = 0, and (if we please) 
A — 1 

We have now to consider the boundary conditions which hold at the 
surfaces of transition. In the case of sound travelling through gas, where 
$ is taken to represent the velocity-potential, these conditions are the 
continuity of clfyjdx and of <j<£, where a is the density. Whether the 
multiplier attaches to the dependent variable itself or to its derivative is of 
no particular significance. For example, if we take a new dependent 
variable i/r, equal to crfa the above conditions are equivalent to the con- 
tinuity of i/r and of a^dtyjdx. Nor should we really gain generality by 
introducing a multiplier in both places. We may therefore for the present 
confine ourselves to the acoustical form, knowing that the results will 
admit of interpretation in numerous other cases. 

At the first transition x = x± the boundary conditions give 

ai (Bi - Ai) = a 2 (B 2 - A 2 ), <n (Bi + Ai) = a 2 (B 2 + A 2 ). (8) 

If we stop here, we have the simple case of the juxtaposition of two 
media both of infinite depth. Supposing B 2 = 0, we get 

5l = ^ /o-i — ch /ai __ cr 2 /cr 1 -'Cot0 2 Jcot0 1 ,gx 

A x cr 2 /<ri + a 2 /ai a 2 /ai + cot 2 /cot 6\ * 

For a further discussion of (9) reference may be made to ' Theory of 

Sound ' (he. cit.). In the case of the simple gases the compressibilities are 

the same, and <n sin 2 0i = a 2 sin 3 6 2 . The general formula (9) then identifies 

itself with FresneFs expression 

tan(fli — 2 ) qqx 

tan (6»i + 2 ) * V ; 

On the other hand, if a 2 = <n, the change being one of compressibility 
only, we find 

8111(02+01) 

FresneFs other expression. 
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In the above it is supposed that a 2 (and 8 2 ) are real. If the wave be 
incident in the more refractive medium and the angle of incidence be too 
great, a 2 becomes imaginary, say —ia 2 . In this case, of course, the reflection 
is total, the modulus of (9) becoming unity. The change of phase incurred 
is given by (9). In accordance with what has been said these results are at 
once available for the corresponding optical problems. 

If there are more than two media, the boundary conditions at x = x 2 

are — 

a 2 {B 2 e ia ^~^-A 2 e^ ia ^-~ x ^} = a 3 (B 3 -A 3 ), (12) 

a 2 {B 2 e ia ^~ x J + A 2 e~ ia ^- x J} = ^(Bg + Ag), (13) 

and so on. For extended calculations it is desirable to write these equations 
in an abbreviated shape. We set 

B2 — A2 = H2, B2 + A2 = K2, etc., (14) 

cosa 2 (x 2 —Xi) = c 1} i sin a 2 (x 2 —xi) = si, etc. (15) 

a z ja 2 = a 2) cr%l cr 2 = fi 2 , etc. ; (16) 
and the series of equations then takes the form 

Hi = otJS.2, Ki = iSiK 2 , (17) 

C1H2+S1K2 = a 2 H3, S1H2 + C1K2 = /82K3, (18) 

C2H3 + S2K3 = a 3 H 4 , S2H3 + C2K3 = /3 3 K 4 , (19) 

and so on. In the reflection problem the special condition is the numerical 
equality of H and K of highest suffix. We may make 

H = -1, K= +1. (20) 

As we have to work backwards from the terms of highest suffix, it is 
convenient to solve algebraically each pair of simple equations. In this 
way, remembering that c 2 —s 2 = 1, we get 

Hi = ai H 2 , Ki = /3iK 2 , (21) 

H 2 = cia 2 H 3 — 5i/3 2 K 3 , K 2 = — «ia 2 H 3 + ci/3 2 K 3? (22) 

H 3 = ^2^31-14—52^31^4, K 3 = — S2a 3 H 4 + c 2 /3 3 K 4 , (23) 

and sq on. In these equations the c's and the /3's are real, and also the 
as, unless there is " total reflection " ; the s J s are pure imaginaries, with the 
same reservation. 

When there are three media, we are to suppose in the problem of reflection 
that H 3 = -1, K 3 = 1. Thus from (21), (22), 

Hi = — «i (Ci«2 + S1&), Ki = /3i (s 1 a 2 + d/8 2 ) ; 

so that — = — + Hl = Cl (ft^ 2 — ai * 2 } + 8l ( (X2 ^ 1 ~" ^ 2 ) (2&\ 

Ai Ki — Hi ci (j3jfi 2 + «i«2) + Si(ot2/3i + oti/32y 
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If there be no "total reflection," the relative intensity of the reflected 
waves is 

Cl 2 (§1§2 — *l*2f — Si 2 ( gggi — cn/3 2 ) 2 ( 2E\ 

where ci 2 = cos 2 a 2 (^2— %i), —Si 2 = sin 2 a 2 (x 2 —%i). (26) 

The reflection will vanish independently of the values of c\ and s\ } i.e., 
whatever may be the thickness of the middle layer, provided 

(3i/3 2 ~ <*i*2 = 0, ct 2 fii ~~ ^1^2 — ; or /3i = ot lf /3 2 = u 2 , 

since these quantities are all positive. Eeference to (9) shows that these 
are the conditions of vanishing reflection at the two surfaces of transition 
considered separately. 

If these conditions be not satisfied, the evanescence of (25) requires 
that either c\ or s\ be zero. The latter case is realised if the inter- 
mediate layer be abolished, and the remaining condition is equivalent to 
0-3/0-1 r= azja ly as was to be expected from (9). We learn now that, if there 
would be no reflection in the absence of an intermediate layer, its intro- 
duction will have no effect provided a 2 (x 2 —xi) be a multiple of ir. An 
obvious example is when the first and third media are similar, as in the 
usual theory of " thin plates." 

On the other hand, if c h or cos a 2 (x 2 — -x{), vanish, the remaining require- 
ment for the evanescence of (25) is that j3 2 j a 2 = fii/ a\. 



In this case 



/3i — oti __ f3 2 ~-a 2 , 



£i + ai /3 2 + oc 2 

;So that by (9) the reflections at the two faces are equal in all respects. 
In general, if the third and first media are similar, (25) reduces to 

{ffi/«i — «i/fti} 2 sin 2 a 2 (x 2 ~-xi) ,^\ 

4 cos 2 a 2 (x 2 —x 1 ) + {ySi/ai -f «i/^i} 2 sin 2 a 2 (x 2 --x 1 ) ' 

which may readily be identified with the expression usually given in terms 
of (9). 

It remains to consider the cases of so-called total reflection. If this 
occurs only at the second surface of transition, ci\, a 2 are real, while a 3 is a 
pure imaginary. Thus ol\ is real, and & 2 is imaginary; ci is real always, 
and si is imaginary as before ; the /3's are always real. Thus, if we 
separate the real and imaginary parts of the numerator and denominator 
of (24), we get 

Ai C 1 ^ 1 ^ 2 + S 1 a 2 ^i + C 1 aici 2 + S 1 a 1 /3 2 > 

of which the modulus is unity. In this case, accordingly, the reflection 



1911.] Waves through a Stratified Medium, etc. 211 

back in the first medium is literally total, whatever may be the thickness 
of the intermediate layer, as was to be expected. 

The separation of real and imaginary parts follows the same rule when 
a 2 is imaginary, as well as a$. For then a x is imaginary, while a 2 , s± are 
real. Thus si« 2 /3i remains real, and Ci<xict 2 , Si&i/3 2 remain imaginary. The 
reflection back in the first medium is total in this case also. 

The only other case requiring consideration occurs when a 2 is imaginary 
and a% real. The reflection is then total if the middle layer be thick enough, 
but if this thickness be reduced, the reflection cannot remain total, as is 
evident if we suppose the thickness to vanish. The ratios a\, a 2 are now 
both imaginary, while s± is real. The separation of real and imaginary 
parts stands as in (24), and the intensity of reflection is still expressed 
by (25). If we take a 2 = —ia 2 , we may write in place of (25), 

(fiil3 2 — ciici 2 ) 2 cosh 2 a 2 (x 2 — #i)— • (*2(3i — a 1 /3 2 ) 2 sinh 2 a/ (x 2 —Xi) ,™ 
(/3i/3 2 + u 1 a 2 f coslrW (x 2 — #i) — (* 2 /3i + <*i/3 2 ) 2 sinh 2 a 2 ' (x 2 —x 1 ) ' 

When x 2 —xi is extremely small, this reduces to 

(#iff2-~*i*2) 2 or (<rs/<ri — aa/(h) 2 

"(ySi/3 2 + *i<% 2 ) 2 ' (crs/cri + as/a^ 2 ' 

in accordance with (9). 

When on the other hand x 2 —X\ exceeds a few wave-lengths, (29) approaches 
unity, as we see from a form, equivalent to (29), viz., 

(Bi 2 — ai 2 ) (/3 2 2 - * 2 2 ) cosh 2 a 2 (x 2 —xi) + (a 2 /3i — <*i/3 2 ) 2 no 

(/3 X 2 - *{*) (/3 2 2 - a 2 2 ) cosh 2 a 2 f (x 2 ^-x 1 ) + (* 2 /3i + *i/3 2 ) 2 ' { } 

It is to be remembered that in (30), <%i 2 , a 2 2 , ol\ql 2 have negative values. 

The form assumed when the third medium is similar to the first may be 
noted. In this case a\a 2 = 1, /3i/3 2 = 1, and we get from (29) 

(A/fti — *i//3i) 2 sinh 2 a 2 (x 2 —xy) ™. 

(/3i/ ol\ — #i//3i) 2 sinh 2 a 2 (x 2 —Xi) — 4c' 

In this case, of course, the reflection vanishes when x 2 —X\ is sufficiently 
reduced. 

Equations (21), etc., may be regarded as constituting the solution of the 
general problem. If there are m media, we suppose H m = — 1, K m = 1, 
thence calculate in order from the pairs of simple equations H m -i, "Km-i ; H w _ 2 , 
K m - 2 , etc., until Hi and Ki are reached ; and then determine the ratio Bi/Ai. 
The procedure would entail no difficulty in any special case numerically 
given ; but the algebraic expression of Hi and Ki in terms of H OT and K m 
soon becomes complicated, unless further simplifying conditions are intro- 
duced. Such simplification may be of two kinds. In the first it is supposed 
that the total thickness between the initial and final media is small 
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relatively to the wave-lengths, so that the phase-changes occurring within 
the layer are of subordinate importance. In the second kind of simplification 
the thicknesses are left arbitrary, but the changes in the character of the 
medium, which occur at each transition, are supposed small. 

The problem of a thin transitional layer has been treated by several 
authors, L. Lorenz,* Yan Ryn,f Drude,J Schott,§ and Maclaurin.[| A full 
account will be found in 'Theory of Light' by the last named. It will 
therefore not be necessary to treat the subject in detail here ; but it may be 
worthwhile to indicate how the results may be derived from our equations. 
For this purpose it is convenient to revert to the original notation so far as. 
to retain a and a. Thus in place of (17), etc., we write 

<2iHi = &2TI2, criKi = o" 2 K2. (32) 

a 2 (ciH 2 + S1K2) = a 3 H 3 , <r 2 (sJBfe + ciK 2 ) = a ^> etc. (33) 

In virtue of the supposition that all the layers are thin, the c's are nearly 
equal to unity and the s's are small. Thus, for a first approximation we 
identify c with 1 and neglect s altogether, so obtaining 

<2iHi = $ 2 H2 = ••• = &mH m , CTiKi = 0*2X2 = ... = CT m K m . (34) 

The relation of Hi, Ki to H OT , K m is the same as if the transition between 

the extreme values took place without intermediate layers, and the law of 

reflection is not disturbed by the presence of these layers, as was to be 

expected. 

For the second approximation we may still identify the c's with unity, 

while the s's are retained as quantities of the first order. Adding together 

the column of equations constituting the first members of (32), (33), etc., we 

find 

aiHi + a 2 SiK 2 + a 3 %K 3 + . . . + flm-iSm-sKm-i = a m~H-m ; (35) 

and in like manner, with substitution of a for a and interchange of K and H, 

eriKi + o- 2 SiH 2 + . . . + o"m-Ai-2H ffl -i = (r m ~K m . (36) 

In the small terms containing s ? s we may substitute the approximate 
values of H and K from (34). For the problem of reflection we suppose 
H m + K m = 0. Hence 

1 ~L- ^" m V* ^^1 

(37) 



Hi_ 


1 -L. ^" m V* ^^1 
(T\Cl m d m " G"2 


Ki 


&m a l -i , a m s? ^Sl 



* 'Pogg. Ann./ 1860, vol. Ill, p. 460. 
t c Wied. Ann., 5 1883, vol. 20, p. 22. 
% < Wied. Ann., 5 1891, vol. 43, p. 126. 
§ < Phil. Trans., 5 1894, vol. 185, p. 823. 
|| ' Roy. Soc. Proc., 5 A, 1905, vol. 76, p. 49. 
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In (37) si = ia 2 (%2—%i), and so on, so that 

^^ = i [«M* , jg^ = < L^ (38 ) 

the integration extending over the layer of transition. 

One conclusion may 'be drawn at once. To this degree of approximation 
the reflection is independent of the order of the strata. It will be noted 
that the sums in (37) are pure imaginaries. In what follows we shall 
suppose that a m is real. 

As the final result for the reflection, we find 

Bi Ki + Hx 



A} iVi~"jLli 



= Re 1 ' 8 , (39) 



where E = «•«/*! -"»/«! , 



tan 8 = 2 



Cirri 



adx-^ 



(40) 
'a?clx 



ci\cr m a m a\ 



amJ °_. (41) 



CCmG'l Ct\<T 



m 



To this order of approximation the intensity of the reflection is unchanged 
by the presence of the intermediate layers, unless, indeed, the circumstances 
are such that (40) is itself small. If a m /ai = a m /a x absolutely, we have 

E = J i^ Ldz-Z™ f^!i^T (42) 

L^mJ Urn) & J 

and 8 = \tt. This case is important in Optics, as representing the reflection 
at the polarising angle from a contaminated surface. 

The two important optical cases : (i) where a is constant, leading (when 
there is no transitional layer) to FresneFs formula (11), and (ii) where 
<rsin 2 is constant, leading to (10), are now easily tretaed as special examples. 
Introducing the refractive index fi, we find after reduction for case (i) 



£/ __ __ 4-7T eos#i 



Ad (fl m 2 '~~ fill 2 ) 



[(jj^-^dx, (43) 



where X 1? fi h relate to the first medium, fi m is the index for the last medium, 
and the integration is over the layer of transition. The application of (43) 
should be noticed when the layer is in effect abolished, either by supposing 
/x, = fjb m or, on the other hand, fi = fii. 

In the second case (42), corresponding to the polarising angle, becomes 

B = . ^4—= fMzWl^Z^ 2 ) dx . (44 ) 

VOL. LXXXVI. — A. Q 
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In general for this case 



cs,, __ ___47TCOS^i 



COS 2 #1 






Xl (ftj-vi 2 ) (cos 2 ^-^sin 2 ft 

(45) 

The second fraction in (45) is equal to the thickness of the layer of 
transition simply, when we suppose fi = //a. 



f 

>i sin 2 ft J 



77 ,, js// ry 47T COS ft Sm 2 ft J /* /,! £\ 

Further, 8 — 8 = — ^ — x — , (46) 

M flm - A*l CQS 2 ^ _ a_ gin 2 ^ 

the difference of phase vanishing, as it ought to do, when ft = pi, or ft m , or, 
again, when ft = 0. 

It should not escape notice that the expressions (10) and (11) have 
different signs when ft and ft are small. This anomaly, as it must appear 
from an optical point of view, should be corrected when we consider the 
significance of 8" — 8'. The origin of it lies in the circumstance that, in our 
application of the boundary conditions, we have, in effect, used different 
vectors as dependent variables to express light of the two polarisations. For 
further explanation reference may be made to former writings, e.g. " On the 
Dynamical Theory of Gratings/'* 

If throughout the range of integration, fi is intermediate between 
the terminal values /^i, fi m , the reflection is of the kind called positive by 
Jamin. The transition may well be of this character when there is no 
contamination. On the other hand, the reflection is negative, if fi has 
throughout a value outside the range between /xi and fi m . It is probable that 
something of this kind occurs when water has a greasy surface. 

The formulas required in Optics, viz. (43), (44), (45), (46), are due, in 
substance, to Lorenz and van Eyn. The more general expressions (41), (42), 
do not seem to have been given. 

There is no particular difficulty in pursuing the approximation from 
(32), etc. At the next stage the second term in the expansion of the c's 
must be retained, while the s's are still sufficiently represented by the first 
terms. The result, analogous to (37), (38), is 



•d 



l—i <r . 

Hi <T\Oj m J o 



* x a 2 dx 



.(T 



.dx + i — I —dx 



&~ &m 



Jo " 



Ki <r m ai 1 _r^ 



rf«2 cx a rd ~> \**) 

crdx ,dx-{-i— \ <rdx 



o°" 



.a 

o °»* J o 

* ' Boy. Soc. Proc.,' A, 1907, vol. 79, p. 413. 
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in which the terminal abscissae of the variable layer are taken to be and d, 
instead of x x and x m ~i. I do not follow out the application to particular 
cases such as<r = constant, or a sin 2 # = constant. For this reference may be 
made to Maclaurin, who, however, uses a different method. 

The second case which allows of a simple approximate expression for the 
reflection arises when all the partial reflections are small. It is then hardly 
necessary to appeal to the general equations : the method usually employed 
in optics suffices. The assumptions are that at each surface of transition the 
incident waves may be taken to be the same as in the first medium, merely 
retarded by the appropriate amount, and that each partial reflection reaches 
the first medium no otherwise modified than by such retardation. This 
amounts to the neglect of waves three times reflected. Thus 

An interesting question suggests itself as to the manner in which the 
transition from one uniform medium to another must be effected in order to 
obviate reflection, and especially as to the least thickness of the layer of 
transition consistent with this result. If there be two transitions only, the 
least thickness of the layer is obtained by supposing in (48) 

ftl — *1 = fi2 — *2 /^gx 

J3i + et! /3 2 + ^ 2 
and 2a 2 (x 2 —xi) = ?r; (50) 

and this conclusion, as we have seen already, is not limited to the case 
of small differences of quality. In its application to perpendicular 
incidence, (50) expresses that the thickness of the layer is one-quarter 
of the wave-length proper to the layer. The two partial reflections are 
equal in magnitude and sign. It is evident that nothing better than this 
can be done so long as the reflections are of one sign, however numerous the 
surfaces of transition may be. 

If we allow the partial reflections to be of different signs, some reduction 
of the necessary thickness is possible. For example, suppose that there are 
two intermediate layers of equal thickness, of which the first is similar to the 
final uniform medium, and the second similar to the initial uniform medium. 
Of the three partial reflections the first and third are similar, but the second 
is of the opposite sign. If three vectors of equal numerical value com- 
pensate one another, they must be at angles of 120°. The necessary 
conditions are satisfied (in the case of perpendicular transmission) if the 
total thickness (2 1) is -J-X, in accordance with 

Q A 
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The total thickness of the layer of transition is thus somewhat reduced, 
but only by a very artificial arrangement, such as would not usually be 
contemplated when a layer of transition is spoken of. If the progress from 
the first to the second uniform quality be always in one direction, reflection 
cannot be obviated unless the layer be at least %\ thick. 

The general formula (48) may be adapted to express the result appropriate 
to continuous variation of the medium. Suppose, for example, that <r is 
constant, making /3 = 1, and corresponding to the continuity of both <p and 
dfyjdx* It is convenient to suppose that the variation commences at x = 0. 
Then (48) may be written 



Ai 



—2i I adx /r-i \ 

3 <>o , (51) 



CLCC —2i\ adx 

— e 

2a 



a at any point x being connected with the angle of propagation by 
the usual relation (3). In the special case of perpendicular propagation, 
a = 27Tfji/\ifii, fi being refractive index and Xi, fii relating to the first 
medium. 

A curious example, theoretically possible even if unrealisable in experi- 
ment, arises when the variable medium is constituted in such a manner 
that the velocity of propagation is everywhere constant, so that there is no 
refraction. Then a is constant, a = 1, and (48) gives 

1^1 1 d(T „-2iax 



Ai J 2<r 



e~ 2lax . (52) 



Some of the questions relating to the propagation of waves in a variable 
medium are more readily treated on the basis of the appropriate differential 
equation. As in (1), we suppose that the waves are plane, and that the 
medium is stratified in plane strata perpendicular to x, and we usually omit 
the exponential factors involving t and y, which may be supposed to run 
through. In the case of perpendicular propagation, y would not appear 
at all. 

Consider the differential equation 

g + ^ = 0, (53) 

m which (unless k 2 can be infinite) it is necessary to suppose that <£ and 
d<p/dx are continuous ; k 2 is a function of x, which must be everywhere 
positive when the transmission is perpendicular, as, for example, in the case 
of a stretched string. When the transmission is oblique to the strata, 

* These would be the conditions appropriate to a stretched string of variable longi- 
tudinal density vibrating transversely, 
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fc 2 may become negative, corresponding to " total reflection," but in most of 
what follows we shall assume that this does not happen. The continuity of 
</) and d$\dx y even though k 2 be discontinuous, appears to limit the applica- 
tion of (53) to certain kinds of waves, although, as a matter of analysis, the 
general differential equation of the second order may always be reduced to 
this form.* 

In the theory of a uniform medium, we may consider stationary waves . or 
progressive waves. The former may be either 

<f) = A cos kqx cospt, or 4> = B sin kqX sin pt ; 

and, if B = + A, the two may be combined, so as to constitute progressive 

waves 

cj> — Acos(pt±Kox). 

Conversely, progressive waves, travelling in opposite directions, may be 
combined so as to constitute stationary waves. When we pass to variable 
media, no ambiguity arises respecting stationary waves; they are such that 
the phase is the same at all points. But is there such a thing as a pro- 
gressive wave ? In the full sense of the phrase there is not. In general, 
if we contemplate the wave forms at two different times, the difference 
between them cannot be represented by a mere shift of position proportional 
to the interval of time which has elapsed. 
The solution of (53) may be taken to be 

4> = A'f(x) + B' x (x), (54) 

where ^r (x), % (%), are real oscillatory functions of x ; A', B', arbitrary 
constants as regards x. If we introduce the time-factor, writing p in place of 
the less familiar c of (1), we may take 

<f> = A cospt . i/r (x) + B sin pt . % (x) ; (55) 
and this may be put into the form 

<£ = H. cob (pt—0), (56) 

where H cos = A\jr (x), H sin = B^ (x), (57) 

or IP = A 2 [f (0)]'+B» [x (%)]\ (58) 

= tan- v ^-xOg) (59) 

A . yjr (x) 

But the expression for <p in (56) cannot be said to represent in general 
a progressive wave. We may illustrate this even from the case of the 
uniform medium where i/r (x) = cos tcx, % (#) = sin /&£. In this case (56) 
becomes 

(f> . = {A 2 cos 2 kx + B 2 sin 2 kx}% cos < pt— tan -1 / — tan kx ) > . (60) 

* Forsyth's ' Differential Equations,' § 59. 
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If B = HhA, reduction ensues to the familiar positive or negative pro- 
gressive wave. But if B be not equal to + A, (55), taking the form 

<f> = |(A+B)cos(l?£— */a£) + i(A— B) cos (pt -\- kx), 

clearly does not represent a progressive wave. The mere possibility of 
reduction to the form (57) proves little, without an examination of the 
character of H and 8. 

It may be of interest to consider for a moment the character of 6 in (60). 
If B/A, or, say, m, is positive, 6 may be identified with km at the quadrants, 
but elsewhere they differ* unless m — 1. Introducing the imaginary 
expressions for tangents, we find 

0= tev+ Msin2« + JM 2 sin4^ + iM 3 sin6/c^+..., (61) 

where M = ? -^=^ . (62) 

ra + 1 v ; 

When /e is constant, one of the solutions of (53) makes <j> proportional to 

e~ iKX . Acting on this suggestion, and following out optical ideas, let us 
assume in general 

= 7 ]e~ i \ adx i (63) 

where the amplitude t) and a are real functions of x, which, for the purpose 
of approximations, may be supposed to vary slowly. Substituting in (53), 
we find 

* ^ 9 + (>c 2 --a 2 ) v --2ia%~{a i v) = 0. (64) 

ax* ax 

For a first approximation, we neglect d 2 r)/dx 2 . Hence 

k = <x, k^k] = C, (65) 

so that <f> = Ctf-*^-*!^, (66) 

or in real form, <£ = Ck~% cos (pt— fa dx). (67) 

If we hold rigorously to the suppositions expressed in (65), the satisfac- 
faction of (64) requires that d 2 r)fda? = 0, or d 2 /c~~%fdx 2 = 0. With omission 
of arbitrary constants affecting merely the origin and the scale of x, this 
makes k 2 = x~\ corresponding to the differential equation 

* 4 S+4> = 0, (68) 



dx 
whose accurate solution is accordingly 

<j> = CxeW-W (69) 

In (69) the imaginary part may be rejected. The solution (69) is, of 
course, easily verified. In all other cases (67) is only an approximation. 
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As an example, the case where k 2 = n 2 /x 2 may be referred to. Here 
^kcIx = nlogx—e, and (67) gives 

<£ = Gx* cos(pt—nlogx + e) (70) 

as an approximate solution. We shall see presently that a slight change 

makes it accurate. 

Eeverting to (64), we recognise that the first and second terms are real, 

while the third is imaginary. The satisfaction of the equation requires 

therefore that 

a* v = C, (71) 

and that tc 2 = CV 4 — ?? ; (72) 

7} dx J 

while (63) becomes <£ = v e"^h""^. (73) 

Let us examine in what cases rj may take the form Dx r . From (72), 

k 2 = C 4 D-^- 4 '-r(r-l)ar 2 . (74) 

If r = 0, k 2 is constant. If r = 1, k 2 = C 4 D~ 4 «^~" 4 , already considered 
in (68). The only other case in which k 2 is a simple power of x occurs when 
r = |-, making 

h? = (C 4 D- 4 + |)^~ 2 = n 2 x~ 2 (say). (75) 

Here rj = Dxi, G 2 rj~ 2 dx = G 2 jD 2 . log a?— e, and the realised form of (73) is 

%> 

<f> = D^cos {p£~ v / (^ 2 — i)l°g^+ 6 }> (76) 

which is the exact form of the solution obtained by approximate methods 
in (70). For a discussion of (76) reference may be made to l Theory of 
Sound/ second edition, § 1486. 

The relation between a and rj in (71) is the expression of the energy 
condition, as appears readily if we consider the application to waves along 
a stretched string. From (53), with restoration of e ipt , 

CLu CLX \^CLX j 

If the common phase factors be omitted, the parts of dfyjdt and d<f>jdx 
which are in the same phase are as prj and a?/, and thus the mean work 
transmitted at any place is as arf. Since there is no accumulation of energy 
between two places, ay 2 must be constant. 

When the changes are gradual enough, a may be identified with k, and 
then i) cc k~%, as represented in (67). 

If we regard t] as a given function of x } a follows when C has been chosen, 
and also k 2 from (72). In the case of perpendicular propagation k 2 cannot be 
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negative, but this is the only restriction. When tj is constant, h? is constant ; 
and thus if we suppose r^ to pass from one constant value to another through 
a finite transitional layer, the transition is also from one uniform k 2 to 
another ; and (73) shows that there is no reflection hack into the . first 
medium. If the terminal values of r\ and therefore of ic 2 he given, and the 
transitional layer be thick enough, it will always be possible, and that in an 
infinite number of ways, to avoid a negative tc 2 , and thus to secure complete 
transmission without reflection back ; but if with given terminal values the 
layer be too much reduced, k 2 must become negative. In this case reflection 
cannot be obviated. 

It may appear at first sight as if this argument proved too much, and that 
there should be no reflection in any case so long as tc 2 is positive throughout. 
But although a constant ij requires a constant k 2 , it does not follow con- 
versely that a constant k 2 requires a constant ?/, and, in fact, this is not true. 
One solution of (72), when tc 2 is constant, certainly is rj 2 = G 2 /k ; but the 
complete solution necessarily includes two arbitrary constants, of which C is 
not one. From (60) it may be anticipated that a solution of (72) may, be 

V 2 = A 2 cos 2 kx + B 2 sin 2 kx = § (A 2 + B 2 ) + J (A 2 - B 2 ) cos 2/ex. (77) 

From this we find on differentiation 

7ff??+*; 2 ^ = * 2 A 2 B 2 ; 

ax* 

and thus (72) is satisfied, provided that 

^a^B 2 = C 2 . (78) 

It appears then that (77) subject to (78) is a solution of (72). The 
second arbitrary constant evidently takes the form of an arbitrary addition 
to x, and rj will not be constant unless A 2 = B 2 . 

On the supposition that t] and a are slowly varying functions, the 
approximations of (65) may be pursued. We find 
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The retardation, as usually reckoned in optics, is ^Kdx. The additional 
retardation according to (80) is 
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As applied to the transition from one uniform medium to another, the 
retardation is less than according to the first approximation by 

(81) 
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The supposition that r\ varies slowly excludes more than a very small 
reflection. 

Equations (79), (80) may be tested on the particular case already referred 
to where tc = n/x. We get 

\ Ion 2 / x\ 8n/ 

;.so that [a dx = I n — — ~ ) log,?;. 

J \ 8nJ c 

When %~ 4 is neglected in comparison with unity, n—^n" 1 may be identified 
with \/(n 2 —\). 

Let us now consider what are the possibilities of avoiding reflection when 
the transition layer (#2— #i) between two uniform media is reduced. If 
'Vi> K i \ Vz> K 2 ar ^ the terminal values, (72) requires that 

We will suppose that r} 2 ^>rji' If the transition from ?;i to ^ 2 be made 
too quickly, viz., in too short a space, d 2 7)/dx 2 will become somewhere so 
large as to render k 2 negative. The same consideration shows that at the 
beginning of the layer of transition (x\), drj/dx must vanish. The quickest 
admissible rise of tj will ensue when the curve of rise is such as to make 
k 2 vanish. When rj attains the second prescribed value 1)2, it must suddenly 
become constant, notwithstanding that this makes /c 2 positively infinite. 

From (72) it appears that the curve of rise thus defined satisfies 

f| = CV 3 . (82) 

dx 2 

The solution of (82), subject to the conditions that tj = 771, drjjdx = 0, 
when x = x±, is 

V 2 - Vl 2 =.CV 2 («-^i) 2 = /eiV(^-^i) 2 . (83) 

Again, when 7) = r} 2 , x = a^, so that 

^(^-^ = 3^ = ^^ a , (84) 

giving the minimum thickness of the layer of transition. 

It will be observed that the minimum thickness of the layer of transition 
necessary to avoid reflection diminishes without limit with #1— -#2, that is, as 
the difference between the two media diminishes. However, the arrange- 
ment under discussion is very artificial. In the case of the string, for 
•example, it is supposed that the density drops suddenly from the first 
uniform value to zero, at which it remains constant for a time. At the end 
•of this it becomes momentarily infinite, before assuming the second uniform 
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value. The infinite longitudinal density at x 2 is equivalent to a finite load 
there attached. In the layer of transition, if so it may be called, the string 
remains straight during the passage of the waves. 

If, as in the more ordinary use of the term, we require the transition to, 
be such that k 2 moves always in one direction from the first terminal value 
to the second, the problem is one already considered. The minimum 
thickness is such that tc 2 has throughout it a constant intermediate value, 
so chosen as to make the reflections equal at the two faces. 

It would be of interest to consider a particular case in which h? varies 
continuously and always in the one direction. As appears at once from (72) r 
d 2 rj/dx 2 , as well as drj/dx, must vanish at both ends of the layer, and there 
must also be a third point of inflection between. If the layer be from x = 
to x = /3, we may take 

^~?= Ax(x-a)(x-l3). (85) 

dx 2. 

We find that /3 = 2 a, and that 



■< ■ 



7j2-Vi 4a 5 120 4 3 J 
_____ P = ~ (X- a) (fc- 2 a). (87) 

rj 2 —y]\dx 2 4a 5 v /v / \ /■ 

From these k 2 would have to be calculated by means of (72), and one 
question would be to find how far a might be reduced without interfering 
with the prescribed character of k 2 . But to discuss this in detail would lead 
us too far. 

If the differences of quality in the variable medium are small, (72) 
simplifies. If t) , k be corresponding values, subject to k 2 = C 4 t;q~ 4 , we- 
may take 

V = VoW> * 2 = *;(> 2 + Stf 2 , (88) 

where tf and S/c 2 are small, and (72) becomes approximately 



dx% 



+ 4zk 2 v = -VoSk 2 - (89) 



Eeplacing x by t, representing time, we see that the problem is the same 
as that of a pendulum upon which displacing forces act ; see ' Theory of 
Sound/ § 66. The analogue of the transition from one uniform medium, 
to another is that of the pendulum initially at rest in the position of 
equilibrium, upon which at a certain time a displacing force acts. The 
force may be variable at first, but ultimately assumes a constant value. If 
there is to be no reflection in the original problem, the force must be of 
such a character that when it becomes constant the pendulum is left at rest 
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in the new position. If the object be to effect the transition between the 
two states in the shortest possible time, but with forces which are restricted 
never to exceed the final value, it is pretty evident that the force must 
immediately assume the maximum admissible value, and retain it for such 
a time that the pendulum, then left free, will just reach the new position 
of equilibrium, after which the force is reimposed. The present solution 
is excluded, if it be required that the force never decrease in value. Under 
this restriction the best we can do is to make the force assume at once half 
its final value, and remain constant for a time equal to one-half of the free 
period. Under this force the pendulum will just swing out to the new 
position of equilibrium, where it is held on arrival by doubling the force. 
These cases have already been considered, but the analogue of the pendulum 
is instructive. 

Kelvin* has shown that the equation of the second order 

£/l^A = yi (90) 

dx\PdxJ u v 7 

can be solved by a machine. It is worth noting that an equation of the 
form (53) is solved at the same time. In fact, if we make 

*>-£• *»-£• (91 > 

we get on elimination either (90) for y ly or 

gf = ^ 3 (92) 

for y 2 . Equations (91) are those which express directly the action of the 
machine. 

It now remains to consider more in detail some cases where total reflection 
occurs. When there is merely a simple transition from one medium (1) to 
another (2), the transmitted wave is 

fa = A 2 e- ia ^ x - x ^e i < ct+ ^\ (93) 

If there is total reflection, a 2 becomes imaginary, say — ia% ; the trans- 
mitted wave is then no longer a wave in the ordinary sense, but there 
remains some disturbance, not conveying energy, and rapidly diminishing 
as we recede from the surface of transition according to the factor e" 0,2 ' ^ x ~ Xl K 

From (2) 

a 4tt 2 2a 4tt 2 /V 1 2 . 2Z) n 

or a 2 ' ^^vW^i-VxVVs 2 ). . (94) 

* ' Koy. Soc. Proc.,' 1876, vol. 24, p. 269. 
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It appears that soon after the critical angle is passed/ the disturbance in 
the second medium extends sensibly to a distance of only a few wave- 
lengths. 

The circumstances of total reflection at a sudden transition are thus very 
simple; but total reflection itself does not require a sudden transition, and 
takes place however gradual the passage may be from the first medium 
to the second, the only condition being that when the second is reached 
the angle of refraction becomes imaginary. From this point of view total 
reflection is more naturally regarded as a sort of refraction, reflection 
proper depending on some degree of abruptness of transition, Phenomena 
of this kind are familiar in Optics under the name of mirage. 

In the province of acoustics the vagaries of fog-signals are naturally 
referred to irregular refraction and reflection in the atmosphere, due to 
temperature or wind differences ; but the difficulty of verifying a suggested 
explanation on these lines is usually serious, owing to our ignorance of the 
state of affairs overhead.* 

The penetration of vibrations into a medium where no regular waves can 
be propagated is a matter of considerable interest ; but, so far as I am 
aware, there is no discussion of such a case, beyond that already sketched, 
relating to a sudden transition between two uniform media. It might have 
been supposed that oblique propagation through a variable medium would 
involve too many difficulties, but we have already had opportunity to see 
that, in reality, obliquity need not add appreciably to the complication of 
the problem. 

To fix ideas, let us suppose that we are dealing with waves in a membrane 
uniformly stretched with tension T, and of superficial density p, which is 
a function of x only. The equation of vibration is (' Theory of Sound/ 
§ 194) 

\dx 2 ' difl p dt* 9 

or, if <£ be proportional to ^'(^+^ as j n ^^ 

^ + V> p /T-V>)4> = 0, (95) 

agreeing with (53) if k? = c 2 p/T-I) 2 . (96) 

* A.n observation during the exceptionally hot weather of last summer recalled my 
attention to this subject. A train passing at high speed at a distance of not more than 
150 yards was almost inaudible. The wheels were in full view, but the situation was 
such that the line of vision passed for most of its length pretty close to the highly heated 
ground. It seemed clear that the sound rays which should have reached the observers 
were deflected upwards over their heads, which were left in a kind of shadow. 
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The waves originally move towards the less dense parts, and total 
reflection will ensue when a place is reached, at and after which tc 2 is- 
negative. The case which best lends itself to analytical treatment is when 
p is a linear function of x. k 2 is then also a linear function ; and, by 
suitable choice of the origin and scale of x, (95) takes the form 

gf +|^ = 0. (97) 

The waves are now supposed to come from the positive side and are totally 
reflected at x = 0, The coefficient and sign of x are chosen so as to suit 
the formulae about to be quoted. 

The solution of (97), appropriate to the present problem, is exactly the 
integral investigated by Airy to express the intensity of light in the 
neighbourhood of a caustic* The line x = is, in fact, a caustic in the 
optical sense, being touched by all the rays. Airy's integral is 

W = cos|7r(^ 3 — miv)dw. (98) 

Jo 

It was shown by Stokesf to satisfy (9*7), if 

x (in his notation n) = {^irf^m. (99) 

Calculating by quadratures and from series proceeding by ascending powers 
of m, Airy tabulated W for values of m lying between m = ± 5*6. For 
larger numerical values of m another method is necessary, for which 
Stokes gave the necessary formulae. Writing 

<fi = 2 QxyP = it (}mfl 2 , (100) 

where the numerical values of m and x are supposed to be taken when 
these quantities are negative, he found when m is positive 

W = 2* (3m)-* {Kcos(£-i7r) + Ssm(£-iir)}, (101) 

1 -|-v -« X.ts.f.XX |X.cr*i*XX«Xc)*Xi.Xt7.^iJc) /"I f\ c\\ 

where a = l-j - 2 (72( ^ )2 + O.S .4(72«£)* " " - ' {WZ) 

q^ X • k) X < tJ • i * XX • X<-J • X I . • ,« /% O \ 

b== 17720 1.2.3(720)3 + **" ilAJ ^ 

When m is negative, so that W is the integral expressed by writing — m 
for m in (98), 

w = a-ics.o^li-f^+f^^-...} • (t04> 

* 'Carnb. Phil. Trans./ 1838, vol. 6, p. 379; 1849, vol. 8, p. 595. 

f ' Camb. Phil. Trans./ 1850, vol. 9 ; ' Math, and Phys. Papers/ vol. 2, p. 328. 

% Here used in another sense. 
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The first form (101) is evidently fluctuating. The roots of W=0 are 
given by 

jli ' n.o^O'O 281 ^ 0*026510, /inK . 

^-•-O^-j^^ (105) 

i being a positive integer, so that for i = 2, 3, 4, etc., we get 

m = 4-3631, 5-8922, 7*2436, 84788, etc. 

For i = 1, Airy's calculation gave m = 24955. 

The complete solution of (97) in series of ascending powers of x is to be 
obtained in the usual way, and the arbitrary constants are readily determined 
by comparison with (98). Lommel* showed that these series are expressible 
by means of the Bessel's functions Ji, J_i. The connection between the 
complete solutions of (97), as expressed by ascending or by descending semi- 
convergent series, is investigated in a second memoir by Stokes.f A 
reproduction of the most important part of Airy's table will be found in 
Mascart's ' Optics ' (vol. 1, p. 397). 

As total reflection requires, the waves in our problem are stationary as 
regards x. The realised solution of (95) may be written 

<£ = Wcoa(ct + by), (106) 

W being the function of x already discussed. On the negative side, 
when x numerically exceeds a moderate value, the disturbance becomes 
insensible. 

*. ' Studien uber die Bessel'schen Functionen,' Leipzig, 1868. 
f < Camb. Phil. Trans.,' 1857, vol. 10, p. 106. 



